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Abstract. When statistical inversion of a lag profile is used 
to determine an incoherent scatter target, the posterior vari- 
ance of the estimated target can be used to determine how 
well a certain set of transmission codes perform. In this work 
we present an incoherent scatter radar transmission code op- 
timization search method suitable for different modulation 
types, including binary phase, polyphase and amplitude mod- 
ulation. We find that the combination of amplitude and phase 
modulation provides better performance than traditional bi- 
nary phase coding, in some cases giving better accuracy than 
alternating codes. 



1 Introduction 

Incoherent scatter radar lag profile measurements can be de- 
convolved using statistical inversion with arbitrary range and 
time resolution as shown by |Virtanen et al. ( 2007a| l. In this 
case, the transmission does not need to be completely free 
of ambiguities. The only important factor is the variance of 
the estimated target autocorrelation function. Even though 
alternating codes are transmission sequences that are opti- 
mal in terms of posterior variance when integrated over the 
code transmission cycle ( |Lehtinen] |1986| l, shorter and only 
slightly less optimal code groups are beneficial in many cases 
where an alternating code sequence is too long. Also, a 
shorter code group offers much more flexibility when de- 
signing radar experiments, e.g., making it easier to combine 
multiple different experiments in the same frequency chan- 
nel and simplifying ground clutter removal. The use of short 
transmission codes is described in more detail in the compan- 



ion paper by Virtanen et al. ( 2007b i submitted to this same 
issue. 



We have previously studied the target estimation variance 
of a coherent target where the target backscatter is assumed 
to stay constant while the transmission travels through the 
target (Vierinen et al. 2006) 1 . We found using an optimiza- 
tion algorithm that a combination of amplitude and arbitrary 
phase modulation can achieve very close to optimum cod- 
ing (the order 10~ 4 less than optimal in terms of normalized 
variance). In this study we apply an optimization method, 
similar to the one used for coherent targets, to find transmis- 
sion codes that minimize the variance of incoherent target 
autocorrelation function estimates. We compare results of 
the optimization algorithm for several different modulation 
methods. 

All formulas in this paper use discrete time, unless oth- 
erwise stated. All waveforms discussed are complex valued 
baseband signals. The ranges will be defined as round-trip 
time for the sake of simplicity. 



2 General transmission code 

A code with length L can be described as an infinite length 
sequence with a finite number of nonzero bauds with phases 
and amplitudes defined by parameters (f>k and a^. These pa- 
rameters obtain values <pk € [0, 2ir] and afc £ [a m i n ,a max ], 
where fc € [1, . . . , L] : k € N. The reason why one might 
want to restrict the amplitudes to some range stems from 
practical constraints in transmission equipment. Usually, the 
maximum peak amplitude is restricted in addition to average 
duty cycle. Also, many systems only allow a small number 
of phases. The commonly used binary phase coding allows 
only two phases: <fik € {0,7r}. 

By first defining 5(t) with tgZas 
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we can describe an arbitrary baseband radar transmission 
envelope e(t) as 



e(t) = ^2a k e i<t,k 8{t-k + l). 



(2) 



fe=i 



We restrict the total transmission code power to be con- 
stant for all codes of similar length. Without any loss of gen- 
erality, we set code power equal to code length 



£K*)I S 



(3) 



t=i 



This will make it possible to compare estimator variances of 
codes with different lengths. It is also possible to compare 
codes of the same length and different transmission powers 
by replacing L with the relative transmission power. 

3 Lag estimation variance 

We will only discuss estimates of the target autocorrelation 
function a T (r) with lags r that are shorter than the length of 
a transmission code (here r is the range in round-trip time, 
and it is discretized by the baud length). The lags are as- 
sumed to be non-zero multiples of the baud length of the 
transmission code. Autocorrelation function estimation vari- 
ance is presented more rigorously in the companion paper 



by Lehtinen et al. (2007 1 submitted to the same special issue. 
The variance presented there also includes pulse-to-pulse and 
fractional lags, taking into account target post-integration as 
well. 

Lag profile inversion is conducted using lagged products 
for the measured receiver voltage, defined for lag t as 



m T (t) = m(t) m(t - 



(4) 



As more than one code is used to perform the measure- 
ment, we index the codes with c as e c (t). For convenience, 
we define a lagged product of the code as 



e c T (t) = e c (t) e^t + r). 



(5) 



With the help of these two definitions, the lagged product 
measurement can be stated as a convolution of the lagged 
product of the transmission with the target autocorrelation 
function 



(6) 



The equation also contains a noise term £ T (t), which is rather 
complicated, as it also includes the unknown target a r (r). 
This term is discussed in detail, e.g., by Huuskonen and 



Lehtinen ( 19961. In the case of low SNR, which is typical 
for incoherent scatter measurements, the thermal noise dom- 
inates and £r(i) can be approximated as a Gaussian white 
noise process defined as 



(7) 



where s 2 is the variance of the measurement noise. 

In this case, the normalized measurement "noise power" 
of lag r can then be approximated in frequency domain as 



Pr 
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(8) 



where e£ (u>) = Tjj {e° (t)} is a zero padded discrete Fourier 
transform of the transmission envelope with transform length 
M ^> L. N c is the number of codes in the transmission 
group and L is the number of bauds in a code. Each code in 
a group is assumed to be the same length, 
will not 

For alternating codes of both Lehtinen (1986) and Sulzer 
( 1993 ) type, P T = 1 for all possible values of r. For constant 



amplitude codes, this is the lower limit. On the other hand, 
if amplitude modulation is used, this is not the lower limit 
anymore, because in some cases more radar power can be 
used on certain lags, even though the average transmission 
power is the same. 

To give an idea of how phase codes perform in general, 
Fig. [T] shows the mean lag noise power for random code 
groups at several different code and code group lengths. It 
is evident that when the code group is short and the code 
length is large, the average behaviour is not very good. On 
the other hand, when there is a sufficient number of codes in 
a group, the performance is fairly good even for randomly 
chosen code groups. Thus, we only need to worry about per- 
formance of short code groups that are sufficiently long. 



4 Code optimization criteria 

Nearly all practical transmission code groups result in such 
a vast search space that there is no possibility for an exhaus- 
tive search. As we cannot yet analytically derive the most 
optimal codes, except in a few selected situations, we must 
resort to numerical means. The problem of finding a trans- 
mission code with minimal estimation variance is an opti- 
mization problem and there exist a number of algorithms for 
approaching this problem numerically. 

A typical approach is to define an optimization criteria 
f(x) with a parameter vector x. The optimization algorithm 
then finds x m i n that minimizes f(x). In the case of transmis- 
sion code groups, x will contain the phase <\> c k and amplitude 
a c k parameters of each code in the code group 



= {a\,. 
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(9) 



There are many different ways to define f(x) in the case of 
transmission code groups, but a trivial one is a weighted sum 
of the normalized lag power P T , with weights w T selected in 
such a way that they reflect the importance of that lag 



(10) 
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In this paper, we set w T = 1 for all lags. This gives each 
lag an equal importance. This is a somewhat arbitrary choice 
of weights, in reality they should be selected in a way the 
reflects the importance of the lag in the experiment. 



5 Optimization algorithm 

As our search method will also have to work with codes that 
have a finite number of phases, we needed an algorithm that 
could also work with situations were an analytic or numer- 
ical derivative of f(x) cannot be defined. We developed an 
algorithm for this specific task, which belongs to the class of 



random local algorithms (Lewis and Papadimitriou 1997 1. 
Another well known algorithm belonging to the same class 
is simulated annealing (Kirkpatric k et al.[ [1983 ), which has 
some similarities to the optimization algorithm that we used. 

The random local optimization algorithm is fairly efficient 
at converging to a minima of f(x) and it can also to some 
extent jump out of local minima. In practice, it is faster to 
restart the optimization search with a different random initial 
parameter set, in order to efficiently locate the minima of 
f(x) that can then be compared the different optimization 
runs. 

A simplified description of our code search algorithm that 
searches for local minima of f(x) is as follows: 

1. Randomize parameters in x . 



1. For a sufficient number of steps, randomize a new value 
for one of the elements of x and accept the change if 
f(x) is improved. 

3. Randomize all parameters x, accept the change if f(x) 
is improved. 

4. If sufficient convergence to a local minima of f(x) has 
been achieved, save x and goto step 1. Otherwise go 
to step 2. The location of the minima can be further 
fine tuned using gradient-based methods, if a gradient is 
defined for f(x). 

In practice, our algorithm also included several tunable 
variables that were used in determining the convergence of 
f(x) to a local minima. Also, the number of local minima to 
search for depends a lot on the number of parameters in the 
problem. In many cases we are sure that the global minima 
was not even found as the number of local minima was so 
vast. 



6 Optimization results 

In order to demonstrate the usefulness of the optimization 
method, we searched for code groups that use three different 
types of modulation: binary phase modulation, polyphase 
modulation, and the combination of amplitude and polyphase 



modulation, which we shall refer to as general modulation. 
In this example, we used a,k = 1 for the constant amplitude 
modulations and allowed amplitudes in the range afe £ [0,2] 
for general modulation codes, while still constraining the to- 
tal transmission code power in both cases to be the same. 

The results are shown in Fig. [T] In this case the re- 
sults are shown in terms of mean lag noise power P = 
(L — J2 T Pt- It is evident that significant improvement 
can be achieved when the code group length is short. For 
longer code groups, the optimized groups do not different 
that much from random code groups. Also, one can see that 
optimized polyphase codes are somewhat better than binary 
phase codes; ultimately general phase codes are better than 
polyphase codes - in some cases the mean lag noise power is 
actually better than unity. The reason for this is that ampli- 
tude modulation allows the use of more power for measuring 
some lags, in addition to allowing more freedom in remov- 
ing range ambiguities. It should also be noted, that when 
the code or code group length is increased, the difference be- 
tween modulation methods also becomes less significant. 



7 Conclusions 

We have introduced an optimization method suitable for 
searching optimal transmission codes when performing lag 
profile inversion. General radar tranmission coding, i.e., 
modulation that allows amplitude and arbitrary phase shifts, 
is shown to perform better than plain binary phase modula- 
tion. Amplitude modulation is shown to be even more ef- 
fective than alternating codes, as the amplitude modulation 
allows the use of more radar power in a subset of the lags. 

For sake of simplicity, we have only dealt with estima- 
tion variances for lags that are non-zero multiples of the 
baud length, with the additional condition that the lags are 
shorter than the transmission pulse length. It is fairly easy to 
extend this same methodology for more complex situations 
that, e.g., take into account target post-integration, fractional 
or pulse-to-pulse lags. This is done by modifying the opti- 
mization criterion f(x). 

In the cases that we investigated, the role of the modulation 
method is important when the code length is short. When 
using longer codes or code groups, the modulation becomes 
less important. Also, the need for optimizing codes becomes 
smaller when the code group length is increased. 

Further investigation of the high SNR case would be ben- 
eficial and the derivation of variance in this case would be 
interesting, albeit maybe not as relevant in the case of inco- 
herent scatter radar. 
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Fig. 1. The mean lag noise power for random codes, optimized binary phase codes, optimized polyphase codes and optimized amplitude 
cifc G [0, 2] and arbitrary phase modulated (general modulation) codes. The largest improvements are achieved for short code groups. Also, 
it is clear that the combination of amplitude and phase modulation provides the best lag variance. 
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